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We explore the properties of the QCD high energy evolution in the limit of a dilute target. Using 
the recently established property of selfduality of the evolution operator (hep-ph/0502119), we show 
how to properly define the target gluon and dipole creation operators in terms of dual Wilson 
lines (dual eikonal factors). We explain how to expand these operators in terms of the functional 
derivatives of the color charge density, in the situation when they act on the eikonal factors of the 
projectile partons. We explicitly calculate the expansion of the high energy evolution operator to 
fourth order in the functional derivatives. Our result is infrared and ultraviolet finite, but does 
not coincide with the formula given in hep-ph/0501088 We resolve this discrepancy by showing 
that the identification of the dipole creation and annihilation operators used in hep-ph/0501088 
is incomplete, and provide the required corrections to these definitions. The use of the corrected 
operators in the calculational framework of hep-ph/0501088 reproduces our result. We also prove 
that there is no discrepancy between the expansion of the JIMWLK equation and the dualization 
of the expansion of the weak field limit. 


I. INTRODUCTION. 

The generalization of the high ener^ QCD evolution to include Pomeron loops has been a subject of vigorous 
research in the last year or so The work Q] lead to realization that the high density limit 

of the QCD high energy evolution equation, the so called JIMWLK equation 113,113,113 does not preserve t-channel 
unitarity and thus the inclusion of Pomeron loops is a must for a fully consistent description of high energy scattering. 
This requires the understanding of the evolution operator of the target wave function not only at high density, but in 
principle at arbitrary density. 

The very interesting series of papers UMM pioneered attempts to derive the low density limit of the evolution, 
the limit opposite to that in which the JIMWLK equation is valid. This goal has been fully achieved in |3. It was 
further shown in Q that the t-channel unitarity together with Lorentz invariance of the scattering matrix require that 
the evolution operator when properly extended to arbitrary density must be selfdual, see below. It was also shown in 
[3 that the duality transformation interchanges the high and low density regimes. This Dense-Dilute Duality (DDD) 
relates the evolution equation derived in |3 to the JIMWLK equation^. 

Despite this impressive progress there is still a number of outstandii^ problems even on the level of comparison of 
the results of different groups. In particular as explained in detail in g the proper account of the noncommutative 
nature of the charge density operators necessitates the introduction of an extra ’’ordering” coordinate, which can be 
identihed with the longitudinal coordinate of the target. The results of |3 and Hi are presented without invoking 
such an extra coordinate. Also the evolution equation of [3 is written in terms of path ordered exponentials of the 
functional derivatives with respect to the color charge density, while the hnal result of [3 (and also of Q) has at most 
four such functional derivatives. In addition, there is some confusion in the literature as to whether the expression 
for the evolution operator for week target with four derivative term can be obtained via the DDD transformation 
from the expansion of the JIMWLK equation to the corresponding order. It has been stated in [3 that this cannot 
be done in apparent contradiction with the proof of [3 of the duality between the full unexpanded expressions. It is 
thus important to establish the correspondence between these different approaches and understand their limitations 
if any. 

The purpose of the present paper is to analyze these questions. Qur main results are the following. We show 
that in the evolution equation derived in [3 the integrations over the longitudinal coordinate x~ can be performed 
explicitly as long as one considers the evolution of the products of the charge density operators which are symmetrized 
with respect to all color indices. Such symmetric combinations arise while averaging the S'-matrix of the projectile 
partons in the wave function of a dilute target. We can thus represent the results of [3 in the same language as [3 
and compare the two expressions directly. The operator of |3 still contains all powers of the functional derivatives, 
as opposed to at most four functional derivatives in |3. The truncation in number of derivatives corresponds to 
the assumption that the projectile as well as the target is dilute (the type of assumption used in the dipole model 
calculations ^3iE3)' The higher derivative terms in the kernel are active if the projectile that impinges on the target 
in question contains a large number of partons which overlap in transverse plane. The expression of |3 is derived 
under assumption that only one projectile dipole can scatter off a given target dipole. From the point of view of the 
DDD of 13, this is a restriction forbidding a target dipole to multiply scatter on the projectile. The more complete 
expression of [3 does not impose such restrictions and allows for arbitrary number of projectile dipoles to scatter on 
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a given dipole in the target, or by duality allows multiple scatteri^s of a given target dipole on a projectile. For a 
meaningful comparison therefore the more complete expression of has to be expanded to the four derivative order. 
Performing this expansion we find that the result of as expected is directly related by the DDD transformation of 
0 to the expansion of the JIMWLK equation. It however does not coincide with the expression given in Q,[3. 

The root of this discrepancy is in the fact that the expressions for the dipole creation and annihilation operators in 
terms of the color charge densities and their derivatives used in are incomplete. The duality transformation allows 
us to define the target dipole (and gluon) creation and annihilation operators exactly. In terms of these operators 
in the large Nc limit, the complete evolution operator of Q indeed has the form of a simple 1 —> 2 dipole creation 
vertex, which is the starting point of the formulation of l^- However expanding these operators in the functional 
derivatives, we find that they have extra terms compared to the operators postulated in Q. The terms absent in 
are precisely those which allow the scattering of two projectile dipoles on a single dipole in the target. Although in 
the kinematics of the dipole model (large but dilute target and projectile) the effect of these terms may be suppressed, 
they should be kept in the evolution if one wants to keep the ” Pomeron splitting vertex” at all away from the dilute 
dipole model limit regime as well as in the regime when the target has a small number of dipoles these 

additional terms contribute to the same order in the evolution equation. Restoring these omitted terms in the result 
of [2 we find complete agreement with the expansion of the evolution operator of |^. We also show that, contrary 
to the statement in there is no discrepancy between the expansion of the evolution equation in the weak field 
limit and the dualized expansion of the JIMWLK equation. The two procedures indeed give the same result if both 
expansions are performed properly. 

This paper is organized as follows. In Section 2, after a brief review of the results of @ and 0, we show how to 
explicitly integrate over the longitudinal coordinate in the evolution equation when the evolution of fully symmetric 
operators is considered. We then perform the expansion of the evolution operator to fourth order in functional 
derivatives. In Section 3 we formulate the dipole model limit of the evolution equation of and show that the 
evolution of the target wave function when expressed in terms of the target dipole degrees of freedom is identical to 
the dipole limit of the JIMWLK equation [i3)l9 viewed from the point of view of the projectile dipoles. We find the 
expansion of the target dipole creation/annihilation operators in terms of the functional derivative of density, and 
identify the missing terms in the expressions of Q . We also comment on the relation between our formalism and that 
advocated in 0. Section 4 contains a short discussion of our results. 


II. THE DENSE-DILUTE DUALITY AND THE EVOLUTION AT LOW DENSITY. 

A. Queen’s croquet party: the setup. 

The evolution equation for the wave function of a small target |T) has been derived in |^. It is conveniently 
formulated in terms of the evolution of the target weight function W defined so that the expectation value of an 
arbitrary operator 0[p(x)\ is given by the functional integral 

{T\d[p{x)\\T) = J Dp‘^W[p{x,x-)]0[p{x,x-)]. (2.1) 

Here p°'{x) is the operator of the color charge density which depends on the gluonic degrees of freedom in the target. 
Note, that although the quantum operator p‘^{x) depends only on transverse coordinate x, the classical variable in the 
functional integral has an additional coordinate x~. This coordinate has to be introduced since the operators p“(x) 
at the same transverse position do not commute. The coordinate x~ can be defined to run from 0 to 1 without any 
loss of generality, and we will use this convention in the present paper. As explained in [^, the weight functional W 
cannot be real as it has to contain a Wess-Zumino term which insures the correct commutation relations of the color 
charge density operators. Its precise form is not important for our purposes. What is important is that it insures 
the correct commutation relations between the charge density operators, and also guarantees that the correlation 
functions of the ’’classical” fields /9“(xj")...p^(x“) do not depend on the exact values of the longitudinal coordinates 
x~, but only on their ordering. We will use this property heavily in the following. 

We will also need the definition of the eikonal factor a(x), which is defined in general through[Il[l3 

a“(x,x")r“ = ^{x-y) {S\y;0,x~) p°‘{y,x~)T°‘ S{y;0,x~)} , S{x;0,x~) = V exp{i J dy~T°‘a°‘{x,y~)} . 

( 2 . 2 ) 

where is the generator of the SU(N) group in the adjoint representation. 
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For a small target this expression can be expanded in powers of p, and to leading order 

a“(x,a;")T“ = ^{x - y) + ... p°‘{y,x~)T°‘, (2.3) 

pX~ nX~ pX~ 

S{x;0,x~) = dy~T^a^{x,y~) — / / dy~dz~T^T^a^{x,y~)a^{x^z~) 

Jo Jo Jy- 

= l + iT'^ [ d?/"a“(x,y“) - i{T“,T^} / dy~a‘^{x, y~) f dz~a''{x,z~) + ... 

Jo ^ Jo Jo 

The physical meaning of the unitary matrix 5'(x; 0,1) is the scattering matrix of a single gluon projectile at transverse 
position X on the target \T). The terms in the expansion of the 5'-matrix en . 112.311 correspond to one- and two gluon 
exchanges of the projectile gluon with the target. In the present paper we restrict our consideration to projectiles 
containing small number of partons. The eikonal factor for such a projectile can be written in the form analogous to 
S{x; 0,1), see 0 


E^p] 


V exp{i f 

Jo 


dy~ 


d^xpl{x)a°-{x,y~ 


)} 


where the projectile color charge density 


(2.4) 


p;{x) = E,T^6ix-x0 (2.5) 

with Xi - the transverse coordinates of the gluons in the projectile. Since the projectile charge density T(x) is a sum 
of delta functions, the effect of the integral over transverse coordinates in eq. iinii is simply to symmetrize the product 
of the 5'-matrices of the individual projectile gluons. 

The total S'-matrix in accordance with eq. EU is given by 

gtotal ^ j W[p{x,x-)] SP[p(x,X-)]. (2.6) 


If the number of the projectile gluons is not too large (parametrically smaller than 0{^) and the target is dilute, 
this expression again can be expanded in powers of a. The exact expression is of no interest to us here, except for an 
obvious property which follows from eq. ESI, namely that it can be written entirely in terms of sums of products of 
the integrated target charge density p(x) = dx~p(x,x~). This will be important in the following. 

The S matrix evolves with rapidity as the target wave function changes under boost. The evolution of the wave 
function in the dilute target limit has been derived in 0. The density correlation functions evolve as 0 


^(p“i(xi)...p“”(x„)) = 
p^(x) [{p^-^ixO+S^iz- 


Ois f {x - z)i {y - 

2 7r Jx,v,z {x - zY iy - 

xi)T“i)...(p“"(x„)+(j2( 


7^(-^{/5“(a;)p“(j/),p“'(a:i)...p“”(x„)} -b 
p^(y)) 


(2.7) 


In this expression the factors p°‘Yxi) are understood as unit matrices in the indices be. 

The meaning of this expression is very simple. The target charge density changes due to the emission of an extra 
gluon in the wave function as a result of boost. The gluon at transverse position z is emitted with probability 
proportional to the square of the Weizsaker-Williams field bf{z) oc f d^x|y 3 yp-p“(x) and contributes an extra adjoint 
charge to the charge density at this position. The second term in eq. E3 is the charge density correlator in the 
component of the evolved wave function which contains this extra gluon, while the first term is the virtual correction 
which restores the normalization. 

The shift of the charge density in eq. (|22I) can be conveniently encoded with the help of the finite shift operator 
. Thus the evolution of the charge density correlators can be written in a more compact form0 


A 

W 


as f {x - z)i {y - z)i 
2 7r Jx,v,z (x - zf {y - zY ^ 

{-\[p^{x), [pYy),r---n) 



WWi 



( 2 . 8 ) 













4 


where the functional derivatives act only on the product of p’s in parenthesis. Translating this evolution equation 
into the evolution equation for the target weight functional W gives directly 0 

AwH = (2.9) 

with the evolution operator 

^ {p“(a:, 0)p“(2/, 0) + p“(a:, l)p“(2/, 1) - 2 p“(a:, 0)i?(z, 0, l)“V(y, 1)} ■ 

( 2 . 10 ) 

where we have defined 


X 


KLWMIJ 


as 

27r^ 




[z - x)i(z - 
(z — x)^(z — 


E(z;0,x-y^ 


-I ab 


"PgJo 


dz-T‘ 


Z-) 


( 2 . 11 ) 


Note, that just like in eq. (EHl) the functional derivatives in R act only on W in eq. (E3 and not on the explicit factors 
p“(0) and p“(l) in Eq. (12.1011 . 

As we will discuss in more detail in the last section, the exponential factor R can be attributed the meaning of the 
eikonal factor for the scattering of the extra gluon of the target on the incoming projectile. It is however important 
to realize that the ordering variable which enters eq. (ITTil has indeed the meaning of x and not x~^ as one might 
naively think. From the point of view of the derivation given above this is quite obvious. The shift operator R shifts 
the charge density of the target by the charge of an extra gluon. As a result the incoming partons scatter also on 
the charge density of this additional gluon. This is most clearly seen if we combine eqs. (ICT and to write the 
evolution equation for the scattering matrix. Integrating by parts and using the fact that y^KLWMij hermitian we 
can write 


^S^otal ^ J jjpa ^KLWMIJ TF[p{x,X-)] . (2.12) 

where now the functional derivatives in R act on the projectile ^'-matrix. Since 

R{z, 0,1) EP[p“(x, x")] = EP[p“(x, x") + r“(52(x - z)] , (2.13) 


we indeed see that the effect of the operator R is to make the projectile scatter off the extra emitted gluon in the 
target. 

We prefer to think of x~ as an ordering variable. However if one does insist on interpreting it as a physical 
longitudinal coordinate in the target wave function, it is also easy to understand qualitatively why the operator R 
is integrated over x~ and not x^. R represents the scattering of the extra gluon of the target emitted during one 
step of the evolution. This gluon is emitted with equal probability at arbitrary longitudinal coordinate x“, and thus 
the eikonal factor has to be integrated over x~ to take into account all possible longitudinal positions of the emitted 
gluon. The other longitudinal coordinate x'^ does not appear in our formulation at all. It has the meaning of the 
ordering variable for the projectile charge density and can be introduced if one wishes to make the projectile charge 
density Pp{x) of eq. 112.511 into a classical commuting variable. We do not have any need for this in the present paper 
and thus do not pursue this route any further. 

Eq. (12.1011 is to be compared with the JIMWLK evolution operator, which describes the evolution of the target 
wave function when the color charge density is parametrically large (p“ = 0{l/as)) 


JIMWLK 


as f {z - x)i (z - y)i 
27r2 .L ,, , {z - x)2 (z - j/)2 
(5 (5 


' x,y,z 

6 


- 2 


6a°-{x,0) 6a°‘{y,0) <5a“(x, 1) <5a“(?/, 1) Sa°-{x,0) 


S{z;0,l) 


ab 


Sa’^iy, 1) 


(2.14) 


The similarity between the two limits of the evolution operator is striking. In fact they are related by a simple 
duality transformation, which we refer to as the Dense-Dilute Duality^ 


Sa ’ 


Sp 


la 


(2.15) 
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It was proved in 0 that the high energy evolution operator valid at arbitrary charge density, which should incorporate 
eg. (12.ion and eg. 112.1411 as its low and high density limits respectively, must be selfdual, that is it has to satisfy 

(2-16) 

The DDD transformation eq. 112.1511 interchanges the high and the low density limits, and this is the origin of duality 
between ^JImwlk ^klwmij^ 

Although expressions for both ^klwmij above involve the ’’ordering” variable x~, we know 

that in the case of JIMWLK this dependence is irrelevant. In particular, if the JIMWLK evolution operator 
^jiMWLK gp|.g fuuctions ot the eikonal factor S'(a:,0,1), it can be rewritten entirely in terms of S and its 

functional derivatives. Let us introduce the left and right rotation generators 

JS(a;) = -tr|5(a:)r“^^|, J£(x) =-tr , Jl{x) = [S{x) Mx)]\ (2.17) 

The relation between left and right rotations follows from the unitarity of matrix S 0 . We can now write 

(2.18) 

and the JIMWLK operator becomes 

^ my) + jr{x) mv) - 2 mx) s{zr mv)} . ( 2 . 19 ) 

Any reference to the ordering variable has disappeared in this equation, it is written entirely in terms of the unitary 
matrix S and functional derivatives with respect to it. 

This leads us to believe that the same should be true for the weak density operator -^klwmij^ rpj^g gj^^ gj? j,gg^ 
of this section is to substantiate this point and rid ourselves of the ordering variable x~. 

In all our calculations from now on we only consider expressions in which the evolution operator acts on functions 
of the integrated charge density 


JIMWLK 


O-s 

2tt'^ 




{z - x)i{z - 
{z — x)^{z — 


F[p{x) = ( dx p{x,x )] (2.20) 

Jo 

As we have noted above, the iS-matrix of any projectile on a dilute target is always a function of this type and therefore 
the set of functions of the type eq. is sufficient for our purposes. 

In the operator language eq. (1?:^ means that we are considering only totally symmetrized correlation functions of 
the charge densities of the type (we omit the transverse coordinates, since the noncommutativity occurs only when 
the transverse coordinates of the charge density operators coincide) 

(r| ^ p“L..p“"|T) (2.21) 

{dl.. .Gn y 

The ordering of the quantum operators in this correlation function obviously does not matter, and therefore the 
relevant correlation function of classical fields will be symmetric under the permutations of the ordering coordinate 
x~. 


(T| ^ p“L..p“"|r) = 

{ai...a„} 


J Dp^W[p{x,x )] Yy )-p“’'(a:„) 


( 2 . 22 ) 


Since the correlation functions of the ’’classical” variables p°'(x~) depend only on the ordering of x“’s and not on 
their values, such correlation function is identically equal to (here our choice of the range of x~ comes in handy, as 
otherwise we would have to normalize the integrals on the right hand side) 


in E 




{dl.. .(In } 


‘|T)= f Dp‘^W[p{x,x )] [ dx^... f )-A“"(a;„) 

J Jo Jo 


(2.23) 
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Thus the correlation functions of the x~ integrals of the ’’classical” field dx~p{x,x~) translate into correlators of 

fully symmetrized combinations of the quantum operators. 

We stress again that this set of functions is sufficient for our present purooses, since scattering amplitude of any 
projectile on a small target involves explicitly only averages of this type^jQ- 

Our aim is first, to understand how to get rid of the ordering coordinate x ~, and second to expand eq. (OUll to fourth 
order in functional derivatives. We have two options to proceed, either in terms of the path ordered exponentials as 
in eg. (12. foil , or directly in terms of the evolution of the correlations of the quantum operators /5“(a;) as in the initial 
setup in . We pursue the first formalism in the body of the paper and give a parallel derivation using the operator 
formalism in the appendix. The results of the two approaches are identical. 


B. Off with her head: getting rid of the ordering variable. 

Our interest is in the averages of the type 

I Dp^ W[p{x-)] F[p] . (2.24) 

where the weight functional W can be written as 

W[p{x,x~)\ = S[i?] a;”)] (2.25) 

The last equation means that the functional Fourier transform of W[p] is a function of the ’’eikonal factor” only, 
namely 




(2.26) 


This must be true for the following reason. As shown in Q, the Fourier transform of the weight function W of a state 
is precisely the 5'-matrix with which this state scatters of some other object. That is to say, if we identify in eq. 
the function a with the eikonal phase for scattering on some hadron H {a being determined by the charge density 
of i?) , then S defined via Fourier transform as in eg. 112.26(1 is precisely the scattering matrix of our target state |T) 
on this hadronic state H. This A-matrix clearly must depend only on the scattering matrix of individual gluons and 
not on any other properties of the eikonal phase a(t)|20|. This is in exact analogy with the statement that in the 
JIMWLK equation we take W to depend only on the matrix S, c.f. discussion around eq. (ITT7I) . Thus eq. 112.2511 is 
merely a statement that the target \T) is a state in th^luonic (or more generally partonic) Hilbert space. A similar 
representation for the weight function W was used in Q (see however discussion in Section 3). 

With the egs. 112.2412. 2 511 we can now proceed to express the kernel terms of p and 5/5p. 

The evolution operator eq. JTTUIl has two structural elements, R{x\ 0,1) and p{x\ 0) (or alternatively p{x] 1)). The 
first one is easily dealt with. Consider the operator R acting on products of 


p{x) 



dx p(x, X ) 


(2.27) 


Obviously 


Sp'^(x, X ) 

Thus it is easy to see that 


and more generally 


^ — p“(y) = c c/ f y ) = - v) = 

dp<^{x, X ) Jo 


R{x; 0, 1) F[p{x)] = exp 


R{x; 0, X ) F[p{x)] = exp 


Sp‘^{x) 


P^{y) 


rj^C 


5p^{x) _ 


x-T<^ 


F[p{x)] 


F[p{x)]. 


6p^(x) 

We conclude that when acting on functions of p the path ordering in R is irrelevant. 


(2.28) 


(2.29) 


(2.30) 
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Now consider the action of p{x] 0) and p{x] 1). It is convenient to think of these operators as acting on the factor 
S[i?] in eg. (12.2511 rather than on the function F. We can then use explicitly the duality between p and 5/5a. 

In full analogy with eqs. ill 712.1 St we have the following relation: 

p“(a:, 0) = 0, 1) p^[x] 1) (2.31) 


Eg. (12.3111 holds as long as both side of it act on an arbitrary function of R. It thus holds inside the integral eq. 

The functional derivatives in R act only on the external factors of charge density (left action), and not on the ones 
which appear in Eq. (loni . This relation as discussed in the appendix can be also verified in the operator formalism 
by explicitly commuting the operator p°^{x) from left to right (or from right to left) through any product of operators 
p°''{xi) in eo. (12.2111 . The sum of multiple commutators that arise in this procedure organizes itself into an exponential 
series which gives eq. (lOTt . _ 

In fact in full analogy with eq. it is clear that p°‘[x\Q) {p°‘{x;x )) is simply the operator of the left (right) 

rotation when acting on R{x; 0, x~), that is 


[p'^(x; 0), R^^(x; 0, a;")] = (T^R(x; 0, x-))^""; [p'^(x; cc"), R^^(x; 0, a;")] = (R(x; 0, x-)T‘^)'’'' (2.32) 

Thus the following more general relation (understood in the same operator sense[^), holds for arbitrary x~ 

p°'{x, 0) = i?“^(a:;0, x~) p^{x, x~) (2.33) 

Utilizing arbitrarines of x~ in eq. llT!Tnil the above equations can be written as 

p^(x, 0) = f dx~ R°‘'’{x]0, x~) p^{x, x~) ; P°'{x, 1) = f dx~ R^ {x\ x~, 1) p^{x, x~) 

Jo Jo 

p°‘{x,x~) = dx/ R°'^{x-,x~,x/) p^{x, x/) (2.34) 

These relations again hold inside the integral eq. _ 

From this point onwards we can use the expressions eq. sequentially to express p(a;; 0) in terms of p{x). We 
start with the expansion 


p°‘{x, 0) = f dxi i?“^(a:;0, x^ ) p^{x, x^ ) = 

Jo 

p“(a;) + [ dx/X/ p^{x, x/) + [ dx/ 

Jo \ J ab Jo 


(2.35) 


2 


RC 


5p‘^{x) 


p\x, x^ ) 


ab 


3 


dx^ 


5p^{x) 


p^{x, x^ ) + 


ah 


The first term in the right hand side is already expressed directly in terms of p{x). Now we ’’unwrap” p{x, Xi ) in the 
rest of the terms using again eq. (ITMIl : 


p“(x, ) = / dx 2 R‘^^{x;x^ , X2 ) p^{x, X2 ) = 

Jo 

p^{x, X 2 ) 


(2.36) 


p“(x) + [ dx 2 {x 2 -xj (t'’ ^ 


5p‘^(x) 


' , - (^2 - f 


dx 


ah 


RC 


5p‘^{x) 


p\x, X 2 ) 


ab 


{^2 -X i 
6 


dX~r, 




p*'(x, X 2 ) + 


After substituting this expression into the right hand side of eq. (12.3511 we find that the term with the first functional 
derivative is expressed entirely in terms of p{x). To find the second derivative term we continue with expansions of 
p{x 2 ) in the rest of the terms and so on. The result can be written in the following form 


p“(x, 0) 


G{T‘= 


5 

5p‘^(x) 


) p\x) 


(2.37) 





















where the function G is defined in terms of the integral of the nested exponentials: 


G{u) 


/ exp{a:i u {exp(x 2 - Xj^)u {expix^ - 2:2 )w...}}} 
Jo 


dx. 


(2.38) 


Although this expression is somewhat implicit, it’s Taylor series expansion is well defined and can be calculated to 
any given order. Up to fourth order the coefficients of the expansion are given by 


Giu) = 


1 + u J dx^ x^ + y dx^ dx 2 + ^1 (^2 ~ ^1 


/ (XX X 

+U^ / dx^ dX 2 dx^ -^- 7^ix2 - a^r) + + Xl{x2 - - X 2 ) 


a;, ^ X ^ 


2 


+ u^ / dx^ dx 2 dxg dx^ + ^(2:2 - a;i ) + ^ x(a:2 - x^ f + {x 2 - x^ )(x 3 - 2:3 ) ^ + 


24 6 

ix3 -X2) + 


1 


(2:2 (2:2 - x - i^y ( x 2 - x -^) _ 2 




2 (2 

(xg -X^y + ix 2 - xC)(x^ - X^)(X4 - X 3 ) 


(2.39) 


Thus to fourth order in derivatives we have (in fact we will only need the following expressions to the third order, but 
we give the last term for completeness) 


p“(x, 0) = 


1 + - T' 


1 


T 


Sp‘^{x) J 12 \ Sp‘^{x) J 720 \ 5p‘^{x) 


1 


rj^C 


p\x) , 


(2.40) 


ab 


or equivalently 


p“(x, 0) = p (x) 


1 - - T' 


4 ( T 


Sp‘^{x) J 12 \ 6p‘^{x) J 720 \ Sp‘^{x) 


1 


rj^C 


(2.41) 


- ba 


The reduction of p(x, 1) is similar. The first conversion is to i?(x; 1, Xg ) rather than i?(x; 0, Xg ). Therefore in the 
first expansion of the analog of eq. we have Xg —> (xg — 1). All the rest of the factors are the same. As a 
result the odd terms in the expansion of p(x, 1) have opposite signs to the corresponding terms in the expansion of 
p(x, 0) while the even terms are the same in both expansions 


p“(x, 1) 




J_ —-— 

720 V Sp^x) 



(2.42) 


Again we stress that all the relations 112.4012.4112.4211 are valid only when appearing in expressions of the form eq. 
with eq. 12.2511 . 


C. Who stole the tarts: the four derivative terms. 


With these expressions we can extract the four derivative term of the evolution operator in order to compare it 
with the result of Q . We rewrite the evolution operator here for convenience 


X 


KLWMIJ 


[p\ =--^ f K{x,y,z) [p“(x, 0)p“(y, 0) + p“(x, l)p“( 2 /, 1) - 2p“(x, 0)Rabiz) p^y, 1)] (2.43) 

Jx.v.z 


First we note, that although the kernel in eq. is the Weiszaker-Wiliams kernel, 

(x - z)i{y - z)i 


K{x,y,z) = 


(x — zy{y — zy 


(2.44) 
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when acting on a wave function which is globally invariant under the color charge it can be substituted by the dipole 
kernel 

M{x, y, z) = - 2^x- zffy - zY 

The argument here is identical to that given in . The kernels K and M differ by terms which do not depend either 
on X or on y. The factors p(a;; 0) and p(x; 1) in ea. l2.4:-{ll are simply the charge density operators acting directly on the 
wave function of the target. Thus adding to K a term that does not depend on x produces an integral over x of the 
charge density operator which acts directly on the wave function of the target. If the target state is invariant under 
the global color charge, this extra term vanishes. A note of caution: the present argument, just like the argument 
of 1^ does not allow one to replace K hy M for all gauge invariant states]^ and therefore in the following we will 
continue using the Weiszacker-Williams kernel K rather than M. However for dipole type states the argument is 
indeed valid and thus for the purpose of comparison with Q the kernels K and M are equivalent. 

We can now use the formulae derived in the previous subsection to find the four derivative terms in eq. (E3S|l- First, 
since there is no third derivative contribution in eq. iionii . no fourth derivative terms are present in the combination 

pYx, 0) p“(y, 0) + 1) pYy, 1) - 2 pYx, 0) p“(y, 1) (2.46) 


In fact this part of the operator contains only two derivatives, and therefore simply reproduces part of the BFKL 
kernel. All higher derivatives therefore reside in the expression 


p°-{x, 0)[i?(z) - YabpHy, 1) = p“(a;) 




1 - I ir 


5p’^{x) 


+ — 1 T' 
12 ' 


5p’^{x) 


5p'^{z) 


rpO 


I - - \T 


Sp^{z) 


i 


+ z: [T 


Sp‘^{z)J 24 \ 5p’^{z) 


12 V 5p<^{y) 


p\y) 


(2.47) 


ah 


Focusing on the four derivative term we find [24 




24 7r2 


K{x,y,z)p^{x) p {y) 


- [Tf 


Tc 


_ rpC 


6pf(x) 


T 


Sp‘^{z) 6p‘^{x) 

S <5 


rj^C 


- T 


5p^{z) 5p‘^{x)\ \ 5p'^{z) 


Sp‘^{z) 

6 


2 r 


_ rjne 


Sp<^{z) Sp’^iy) 

e S _ j.e ^ 

Sp^iz) Sp‘^{y) 


TC 


_ rpC 


Sp‘^{z) 6p‘^{x) 


rj^C 


Sp‘^{z) 


TC 


_ rj^e 


Sp<^(z) Sp^{y) 


Tf 


dpf{y) 


(2.48) 


ab 


This operator possesses the required properties of the ultraviolet and infrared finiteness. The ultraviolet finiteness 
is the fact that if the coordinate z becomes equal to x or y, the integrand in eq. vanishes. The infrared finiteness 
is guaranteed since each individual term in the expression contains at least one derivative with respect to p{z). Thus 
if in the initial wave function there are no color charges beyond some radius L, this radius automatically takes on the 
role of the infrared cutoff in eq. EH- 

We can now compare our result with the expression for the four derivative term given in 


Xmsw — 


16 7r2 m 


K{x,y,z) p“(x) p^{y) 


'x,y,z 


Sp^{z) S p^{x)\ \_5p‘^{z) Sp‘^{y) 


(2.49) 


The two look quite different. To make the comparison more explicit we have to take the large W limit of eq. iinHli . 
This procedure is slightly involved and is discussed in more detail in Appendix B. However to see that the two 
expression indeed do not coincide, we can for example act with the two kernels on four factors of p representing four 
projectile partons. At large Nc the action of on a product of two pairs p°'{xi) p°^{x 2 ) p^ixs) p^{x 4 ) is identical to 
the action of the following kernel 


XNc — 


CXc 


48 7r2 m 


K{x,y,z)p°’{x) p°-{y) 


(2.50) 
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X |2 

-f 

■ J 6 

5 

'6 S 6 ' 

■ J 6 

5p^{z) 5 p^{x) 

S 6 ' 

Sp^{z) 
6 ■ 

5p^{z) Sp'^{x) dp'^{y) 

d ^6 „ <5 ■ 

0 o 

_5p‘^{z) Sp‘={y) 

■ (5 6 

5p>{z) 5 p'^iy) 

Sp^{z) 

5p‘^{z) Sp‘^{x) dp‘^{y) 

dp‘^{z) Sp’^ix) 


Sp^{z) 6p^{x) 


Sp^(z) Sp'^{y)\ Sp‘^{z) 


Sp‘^{z) 5p'^{x) Sp‘^{y) 


This clearly does not coinside with eq. lEnnii- In the next section we explain the origin of the discrepancy. 

We note also that the expansion of our operator eq. iinsii yields also three derivative terms, which are not formally 
subleading in the large Nc limit, the odderon type structure similar to the ones discussed in Such terms are 

also absent in UM We will however concentrate on the four derivative terms for definiteness. 


III. DIPOLES OF THE TARGET. 

The Dense-Dilute Duality between the evolution operator eq. and the JIMWLK operator eq. iimii allows us 

to discuss the dipole limit in a way completely analogous to the discussion of the dipole limit in @ . 

We remind the reader, that in the framework of the JIMWLK equation one can define a ’’projectile dipole creation 
operator” as 

s{x,y) = ■^i^[SFix)Sl,{y)] (3.1) 

where Sp is the fundamental representation of the single gluon scattering matrix S{x; 0,1) defined in eg. ( 12 . 211 . The 
reason one can think of it as a dipole creation operator, is that the presence of an extra dipole in the projectile wave 
function \P) results in an extra power of s in the expression for the scattering matrix averaged over the projectile 
wave function 


^P[S] = {P\S\P) 


(3.2) 


By the same token, applying the functional derivative with respect to s{x, y) to S, diminishes the power of s, and 
thus is identified with the projectile dipole annihilation operator. 

It is then very clear how to define creation/annihilation operators for dipoles of the target. In the original derivation 
of [3, the path ordered exponential R{z) appeared directly as a consequence of the presence of an extra emitted gluon 
in the target wave function at transverse coordinate z. Also as we have discussed earlier (see eas. (l2.25l2.2till l if the 
target is scattered on another hadron, each factor of R turns into an eikonal factor of a target gluon. Thus R{x] 0,1) 
should be understood as the creation operator of a target gluon. The target dipole then is created by 


S{x,y) =r{x, y) = tr[RF{x) R\,{y)] 
Conversely, the dipole annihilation operator is simply 


(3.3) 


d{x,y) 


5 

Sr{x,y) 


(3.4) 


To derive the target dipole model we assume in full analogy with 0 that the target wave function depends on 
dipole degrees of freedom only. Referring to our earlier discussion this means that in eq. we take E = S[r]. In 

this case, we can repeat verbatim the derivation of the dipole limit given in [^. We recall that in the case of JIMWLK 
the squares of the left and right rotation operators to leading order in l/Nc are given by 


jR{x)J'h{y) = Jl{x)Jl{y) = ^ J d^u(fv[S{u-y)-S{v-y)][S{u-x)-S{v-x)] s{u, v) 

JL{x)JR{y) = ^ [ d^u(fvtiiS'<{u)T^S{v)T^] [J(u -y)- S{u - y)] [J(u - x) - 5{u - x)] ^ (3.5) 

Nc J ds[u,v) 

An identical calculation gives in the present case 

p°-{x, 0)p“(?/, 0) = p“(x, l)p“( 2 /, 1) = ^ / d^ucfv [J(m -y)- S{v - y)] [J(m - x) - S(v - x)] r{u, v) —- (3.6) 

2 J dr[u, V) 
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The evolution operator in the dipole approximation becomes 

{x - yf 




as 

2 TT 


'x,y,z 


{x- zY (z - yY 


[-r(a;, y) + r{x,z)r{y,z) ] 


6r{x,y) 


(3.7) 


The first lesson from this calculation is that in terms of the target dipole creation and annihilation operators the 
evolution is exactly the same as postulated in Our evolution operator is simply a (1 ^ 2) splitting vertex for 
the target dipoles. The second thing to note however, is that the expressions of the dipole creation and annihilation 
operators used in are incomplete. In particular, the expressions used in Q are 


^MSW 


{x,y) = 1 + 


4.N, 


6p°'{x) Sp‘^{y) 


- NcdMSw{x,y) = p°‘{x)if{y). 


(3.8) 


On the other hand our creation and annihilation operators are infinite series in powers of functional derivatives of 
p. Clearly, any correction to d) or d of eq. (E3 to order {5/5pY and p^{5/5pY respectively will affect the form of the 
evolution operator eq. It is easy to see that such corrections to both, the creation and annihilation operators 

exist. In particular for the dipole creation operator eq. (lOl yields directly 


d\x,y) = 1 + 


[ (5 <5 ] 

^ 1 

[5 <5 ] 

2 

[ (5 (5 ] 

_Sp°-{x) Sp^{y) 

32A^2 

5p°‘{x) 5p°‘{y)_ 


Sp^ix) Sp’^iy) 


I 


96iV2 
4 


dp'^(x) J \6p‘^{x) 
2 


Sp°-{x) J \Sp‘^{x) Sp‘^{y) 


Sp‘^{y)J \Sp^iy) 
2 


Sp<^{x) 5p<^{y) 


-4 


dp°‘iy)J \Sp‘^ix) 6p'^{y) 




(3.9) 


We do not write out explicitly the three derivative terms in df as they do not enter To find the expression of the 
dipole annihilation operator we compare eq. (imHi with eq. (ICT . On one hand from eq. (ICT we have (as we are only 
interested in the second derivative terms in d, we do not have to keep higher order terms in the expression for rfl in 
the equations below) 


p“(x;0)p“(a;;0) = - Ncr{x,y) 


6r{x,y) 


= -nAi + 


4Nr 


5p°‘{x) dp°‘{y) 


d{x,y) , 


(3.10) 


and on the other hand from eq. innii 


p“(x;0)p“(y;0) = il-i(T' 


5p^{x) 


rpa 


6p°-{y)J 12 


r 


Sp‘^{x) 


+ T 


SpYy) 


p‘'ix)p{y)- (3.11) 


ab 


On balance we find the following relation 


- Nc d{x, y) = < 1 - 


- -\T 


4Nc [Sp‘^{x) Sp‘^{y) 


rpa 


5p^{x) 


1 


Sp°-{y)J 12 


(3.12) 


Tc 


5p^{x) 


TC 


SpYy) 


P‘'{x)p\y). 


- y ab 


We thus clearly see that the relation between the dipole creation/annihilation operators and the color charge density 
assumed in 01 is incomplete. The corrections of 0{S^/SpY to dYx,y) and 0{p^Y jdp^') to d{x,y) not included in the 
expressions of 0,i change the form of the splitting vertex when written in terms of the color charge densities and 
its functional derivatives. 

As we discuss in the next section, the extra terms in d and d^ not kept in 0 are responsible for processe when 
more than one dipoles (gluons) of the projectile scatter on the same dipole (gluon) of the target. In the regime 
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where both the target and the projectile are dilute but large (contain large number of dipoles/gluons but low density 
thereof) the effects of the extra terms we find on the total scattering amplitude are presumably suppressed. However 
if the ” Pomeron splitting” terms are important in the regime where the ” Pomeron merging”, or JIMWLK evolution 
becomes effective, all four derivative terms should be considered on par. 

Finally, it is interesting to establish connection between our approach and that of Ref. 0. In the latter paper the 
evolution operator has been written in terms of the projectile dipole degrees of freedom s and S/Ss, rather than the 
target charge density degrees of freedom p and 5/5p. For a dilute target it is easy to relate the two: 


Thus 


s{x, y) = I - {a^{x) - a‘^{y)f 


(3.13) 


P“(^)P“(J/) 


Nc 

al 


dldls{x,y) 


(3.14) 


When acting on F'js] the 
respect to s: 

5p°‘{x) 5p°-{y) 


second derivative with respect to p in the large Nc limit becomes a single derivative with 

_ f — — — — (3.15) 

2 J ds{u, V) 


Similarly, the four derivative term is converted into two derivatives with respect to s. When written in terms of s the 
four derivative term of 01 becomes 


,,4d _ 2as7r 

XmSW — 2 

^.<5 


K(x, y, z) -f(x,z; u, v) -f{y, z; p, r) dl dl s{x, y) 


5^ 


> x^y,z,u,v,r,p 


Ss{u^ v) 6s{p, r) 


(3.16) 


with 7 standing for the dipole-dipole elastic amplitude and defined through 


j{x,z-,u,v) 


32 7r2 


In^ 



u)^ {z — u)^ \ 

vY {z — uY J 


(3.17) 


After integrating the transverse Laplacians by parts we find that this is equivalent to the 2 —> 1 merging term of 0 
(see Eqs. (2.14) and (3.30) of that paper). This is not surprising since the underlying physical assumptions are the 
same in both papers. We thus conclude that the approach of 0 is limited by the same kinematical range (large but 
dilute projectile/target) as that of Hi 

It was also argued in 0 that the merging term s must be accompanied by the virtual contribution of the 

form s s . This term was introduced to preserve the probabilistic interpretation of the dipole evolution of the 
projectile wave function([^ 1^1 and is necessary to ensure the s-channel unitarity: to be consistent with unitarity 
any dipole evolution kernel x[s] has to vanish both at s = 0 and s = 1. We cannot rule out or confirm the presence of 
such a term within our present framework, since the present derivation explicitly assumes dilute target. In the dilute 
target limit the square of the dipole amplitude is linearly related, and thus indistinguishable from the amplitude itself: 
s(a;, y)s(u, v) = s(x, y) + s{u, v) - I. 


IV. DISCUSSION 

In this paper we continued the study of high energy evolution in the dilute target limit. We have shown how to 
integrate over the ” ordering variable” in the evolution equation in this situation. This allowed us to compare the 
results of 0 with those of 0,0,0. We found that although the expression of 0 in the dipole model limit reduced 
precisely to the usual dipole creation vertex, when translated into the language of functional derivatives with respect 
to p, our result contains extra terms relative to those in 0 , 0 . The origin of these extra terms is the contribution to 
the target dipole creation/annihilation operators due to scattering of two projectile dipoles on a single dipole of the 
target. 

This last point becomes obvious once one realizes that the meaning of the dual path ordered exponent R{z) is the 
eikonal phase for scattering of the complete projectile wave function on a single gluon of the target at transverse 
position z (or equivalently, the eikonal phase of the scattering of one gluon of the target on the projectile). We 
have already seen this in the discussion in section 2. To make this more explicit, consider for example a single-gluon 
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projectile impinging on a single-gluon target. The 5'-matrix of one projectile gluon is given by (we neglect the color 
structure throughout this qualitative discussion) 

^(a;) = 1-h ZQ!{a;) (4.1) 

The one gluon target is described by the ’’weight functional” 

W[p\ = R{z)5{p) (4.2) 

This expression for the weight function is analogous to the one used in MM for the dipole model. The vacuum is 
described by a delta function, as an iS-matrix of any projectile on it is equal to unity. The 5'-matrix for this process 
is given by 

J dpSW = 1 -f (4.3) 

which is precisely the gluon-gluon scattering amplitude due to a single gluon exchange. In this toy calculation only the 
first derivative term in R is active since higher derivatives vanish when acting on the single gluon S matrix eq. EH. 
For the projectile which consists of two gluons, the S'-matrix is 

S{x,y) = (1 -h ia(a;))(l -f ia{y)) (4.4) 

The scattering amplitude on a single gluon of a target is then 

/ I 1 ^11 

dpS{x,y)W = 1 -h i{g'^-^{x,z) + g'^-^{y,z)) - z)^{y, z) (4.5) 

The quadratic term comes from the action of the second derivative term in R on the quadratic term in a in eq. E3|). 
Clearly when we repeat this gedankencalculation for the S matrix of n incoming gluons we find that the n-derivative 
term in the expansion of R becomes active, and gives the amplitude for simultaneous scattering of all n incoming 
gluons on the single gluon of the target. It is straightforward to repeat this argument for the target dipole operator 
r. The four derivative term in the expansion of r{x, y) contributes to the amplitude of the simultaneous scattering of 
two projectile dipoles on the single dmole of the target. These are precisely the terms that distinguish our expression 
for d){x,y) from the one used in 

Finally we note that the derivation presented in this paper answers the query raised in 0 about the apparent 
difference between the dual form of the expansion of derivative term in the weak field 

limit of the evolution. When expanding the JIMWLK evolution operator one has not only to handle carefully 
the dependence of the eikonal factor S on the ordering variable x~, but also relate S/Sa(x,l) and S/6a{x,0) to 
S/S J dx~a{x,x~). Expanding naively only the exponental factors S{x), etc. in the JIMWLK kernel gives a result 
which depends on the exact form of the kernel used (i.e. whether it is expressed in terms of S/Sa{x, 1) or S/Sa{x,0) 
or in a mixed representation). The simple reason is that the remaining factors {S/Sa{x, 1) or S/Sa{x,0)) themselves 
contain corrections in powers of a when expressed in terms of S/S f dx~a{x,x~). The proper way to perform the 
expansion is by employing the dual variant of the procedure described in Section 2 of the present paper. This in 
particular generates a term involving the product of functional derivatives in three different spatial positions. The 
result of the correct expansion is the exact dual of our expressions eas. (12.4812.5011 . 


APPENDIX A: OPERATOR ORDERING APPROACH 


In this appendix we perform the calculations of Section 2 in the language of quantum operator ordering. Let us 
denote a fully symmetrized (with respect to color indices) product of n charge density operators by ■ ■ ■ p®"-}. 
Once the product is fully symmetric, the ordering of the operators is irrelevant. This particular set of correlation 
functions has therefore a representation in terms of a weight function IT[p(x)] of a classical field p{x), where the field 
p{x) does not have a dependence on an additional ordering variable. 


nl 




W[p{x)] 


(AI) 


The rapidity evolution of this correlation function follows directly from eq. (ITKll which is valid for any product of 
the charge density operators, including fully symmmetric combinations 


A 

W 


({/5“ 


P“"}) = 


2 7T 


jx - z)i {x - z)i 
{x — z)2 (x — zY 


(A2) 


[P“(r/), 


’}]]) + (P“(^) 


— 1 


ab 


{r ■ ■ ■ p\y)) 
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The correlator in the right hand side of eq. is not yet given in terms of a fully symmetrized products: the factors 
p°‘{x) and p^{y) are not symmetrized with other operators. Our goal is to reformulate (EH) as an equation for the 
distribution W. In order to achieve this we need to rewrite EH) in terms of fully symmetrized combinations of 
This can be done using the canonical commutation relations of operators p: 

[p%x),p\y)] = T:,p‘={x)S{x - y) (A3) 


In order to illustrate the operator approach let us first rederive the relation ESBl- Consider a string of n operators 
multiplied from the left by the operator p^{y) 

p''{y) 


When such a correlation function is represented in terms of a classical field p{x,x~), as in Section 2, the operator 
p^{y) becomes p^{y, x~ = 0) as its position is to the left of all the other operators in the correlation function. In order 
to obtain p^{y, x~ = 1) we need to commute the operator p^{y) to the right of the string. To see how the exponential 
operator arises in this procedure we first rewrite the commutation relation (IA3II in the following form 


p\y)p‘^{x) = p'^{x)p\y) + p‘'{jj)T^i, 


5p^{y) 


p%x) = p“(cr) 


I + T 


Sp‘^{y)_ 


nv) 


(A4) 


be 


In the last equality in IIA4II . the derivative operator acts to the left. Applying the procedure twice we obtain: 


p\y) r{xi)p^-{x2) = p^-{xi)r{^2) 


I + T“ 


5p<^{y) 


+ ^ [T 


6p‘^{y) 


p^iy) 


(AS) 


be 


After repeatedly applying the procedure n times, the operator p^{y) is moved to the very right of the string: 


pHv) P°'^{xi) •••p“"(a:„) = /5“i(xi) •••p“"(x„) 





This can be written as the following relation 


(A6) 


pliv) = 




6 

6p^{y) 


be 


PRiy) 


(A7) 


which should be understood in the precise sense defined by ESJ. Once the correlation function is represented in 
terms of classical fields p{x,x~), Eq. EH) becomes the relation ESU- Note that by construction, the functional 
derivatives in the exponential factor do not act on its “parent” p'yiy). 

We now return to eq. EH) with the aim to rewrite the right hand side of the evolution equation explicitly in terms of 
the correlation functions of fully symmetrized combinations of /5“”s. The symmetrization procedure is only nontrivial 
when the charge density operators are at the same transverse position. Otherwise the operators commute and can 
be moved through each other without obstruction. We will therefore in the following omit the transverse coordinate 
and assume that all the operators are at the point x or y, depending on the coordinate of the external operator which 
comes from the evolution kernel. 

We begin with the double commutator term in eq. ESI. Using the property: 


[A,BCD...] = [A, B]C D ... + B [A, C] D ... + B C [A, B] ■ ■ ■ + ■ ■ ■ 


we find 


[/5“(a:),{p“^ ■••/)“"}] = ■••(z)...p“"}, (A8) 

i 

where {i) denotes the absence of the i-th term. The right hand side of 1A8II is now written in terms of fully symmetric 
products. We can write this as 

([p“(x), =Y.TaaA{p'^r •••(*)••• /5“"}) = T^aAip'^lr (A9) 

i ^ 

Using the same procedure twice we find for the double commutator 

([p“(:r), [pAy)Ar---AA]]) = TATt{{pAAj^^pAy)j^^p^^ (Aio) 
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Thus the double commutator term gives rise to two derivatives with respect to p and no higher derivative terms, for 
any correlator independently of the number of inserted operators n. These are the two derivatives which appear in 
the BFKL limit of Eq. 1A2II . 

We turn now to the symmetrization of the last term in eq. (ESI). We are seeking the result in the form of an 
evolution operator, which can be formulated as a polynomial in derivatives S/6p which act on a symmetrized product 
of operators Our aim is to calculate the terms up to fourth derivative. Since at least one derivative S/dp{z) comes 
from the expansion of the exponent in eq. EH), we only need to extract terms up to the third derivatives that arise 
due to the symmetrization procedure. To extract the terms with up to three derivatives it is sufficient to consider the 
evolution of the one, two and three point functions. 

We consider n = 1, n = 2 and n = 3 in turn. 

• n = 1: 


P P 


= o Pt} + t: 


(All) 


n = 2: 


p‘^\{rr} = \{p^rr} + I 




+ (A12) 


n = 3: 


/3“ r r ^ 


6 

1 

'l2 


P^ P‘^^} 




2 ^ J Lr r J ' 2 ^ J 2 

We immediately infer from ens. ljAl 1IA13ll that up to (and including) the third derivative term we have 


(p“(x){p“^.-.p“"}) = 


1 + ^ ( T’ 


5p^{x) 


+ — ( T 
12 


and similarly 


{{r ■ ■ ■ P^A p\v)) = 


1-T 


5p“(y) 


1 

12 


rpa 


5p^{x) 


5p<^{y) 


{{p^{y)p^^---pA) 


(A13) 


(A14) 


{{P^{y)p^^---P^-)) (A15) 


This coincides with the results obtained in Section 2 with the help of reordering of the path ordered exponentials. 

We have also checked that the operator reordering calculation reproduces the coefficient of the four derivative term 
given in Section 2. The calculation of the four derivative term proceeds along the same lines as above, but the algebra 
is significantly lengthier. As it does not add new understanding we are not reproducing it here in any detail. This 
concludes the explicit demonstration that the operator symmetrization method is equivalent to reordering of path 
ordered exponentials. 


APPENDIX B: LARGE LIMIT. 

In this appendix we discuss the large Nc limit of the kernel -^klwmij^ This is peripheral to the subject of this 
paper, but may nevertheless turn out to be useful in future. As before we concentrate on the four derivative term, but 
the bulk of the discussion applies to any term in the derivative expansion of the kernel. For convenience we rewrite 
the four derivative term 
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Since the kernel does not appear on its own but rather in the integral 

J Dp^ W[p] Fm ■ (B2) 

we should be considering the large Nc limit of this expression. For simplicity we assume that in the expansion of F 
the charge densities are always pairwise contracted in the colour indices, that the only terms present are of the type 


[p^%X2W^y2 )]... (B3) 

This assumption is equivalent to the two gluon exchange approximation for the scattering of each individual dipole 
of the projectile. 

The simplest F of this form which is not annihilated by the four derivatives in y is 


It can be checked by explicit calculation that the action of x'^‘^ on this product of the two pairs is equivalent to 
the action of x'n^ of eq. 11301). However some thought makes it evident that this expression does not produce all 
the leading large Nf, contributions in eq. for more general F. The functional S itself also contains functional 
derivatives. In the simplest term where the number of factors of p in F’ is the same as the number of the functional 
derivatives in any functional derivative from S is forced to act on p in y. For all the terms of this type eq. 
gives the leading large contribution. However if we consider a term in F containing for example a product of three 
pairs of p’s, two of the p factors will be annihilated by the derivatives coming from S. The colour structure in these 
terms will be different and the action of eq. 12.5011 does not emulate correct large Nc limit for these terms. 

We did not find a more simple and compact expression for y^^ which would reproduce the large Nc limit for an 
arbitrary term in F, than simply follow the diagrammatic rules by’t Hooft. This amounts to the following. The group 
SU{N) is promoted to U{N). Every adjoint index a is promoted to a pair of a fundamental and an antifundamental 
indices a ^ (g), and every generator Tj“ is represented by an antisymmetric combination of two vertices 


rpa _ 


1 

V2 






(B5) 


The charge density operator p also acquires a pair of fundamental/antifundamental indices p^ ^ p%. Diagramatically 
then the kernel can be represented as a sum of the terms of the type (we draw a diagram that corresponds to an 
arbitrary set of four transverse coordinates in the functional derivatives - in the actual kernel those coordinates become 
X, y OT z depending on which term in eq. one is representing) To see how this produces the leading large Nc 


FIG. 1: Double line representation of different terms in the kernel y^'’*. 

result for ea. llB2ll we also have to represent diagramatically F and W. The contracted pairs of the charge densities are 
simply represented (Fig. as double wedges (the angle of the wedge is where the indices of the two charge densities 
are contracted) Finally to discuss W we take it in the form eg. (12.2611 with S having a similar structure to F, that is 
where functional derivates come in pairs with pairwise contracted color indices. Diagrammatically each such pair in 
W are therefore represented as in Fig. |21b. 

With these graphical elements one can calculate the large Nc limit for any particular term in F and W. For 
example taking two pairs of p in E and one pair of functional derivatives in W we have the diagrams (Fig. |2I) which 
reproduce the action of eq. E3nil . Other terms can be diagramatically represented, and as usual the leading large Nc 
contribution comes from planar diagrams. To illustrate the fact noted in Section 2 that eq. does not resumm 

all the leading 1/Nc terms we give and explicit diagrammatic example in (Fig. In Fig. 2|the kernel y'*'^ is taken 
to act on six factors of p (3 pairs) out of which four are annihilated by the kernel itself while the remaining ones are 
removed by the action of two derivatives in E of target. Diagramms of this type are not contained in eq. E3ni) . 
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FIG. 2: Double wedge representation of (a) the projectile charge 


density, (b) the target charge density. 



FIG. 3: The planar diagrams contributing to the large Nc limit of the quartic term in F. These diagramms are resummed in 
ea. l2.5(lil . 
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